THE INITIAL-BOUNDARY VALUE PROBLEM FOR THE 
COMPRESSIBLE VISCOELASTIC FLUIDS 



XIANPENG HU AND DEHUA WANG 

Abstract. The global existence of strong solution to the initial-boundary value problem 
of the three-dimensional compressible viscoelastic fluids near equilibrium is established in 
a bounded domain. Uniform estimates in W^''' with q > 3 on the density and deformation 
gradient are also obtained. All the results apply to the two-dimensional case. 



1. Introduction 

Elastic solids and viscous fluids are two extremes of material behavior. Viscoelastic 
fluids show intermediate behavior with some remarkable phenomena due to their "elastic" 
nature. These fluids exhibit a combination of both fluid and solid characteristics, keep 
memory of their past deformations, and their behaviour is a function of these old defor- 
mations. Viscoelastic fluids have a wide range of applications and hence have received 
a great deal of interest. Examples and applications of viscoelastic fluids include from 
oil, liquid polymers, mucus, liquid soap, toothpaste, clay, ceramics, gels, some types of 
suspensions, to bioactive fluids, coatings and drug delivery systems for controlled drug 
release, scaffolds for tissue engineering, and viscoelastic blood flow past valves; see [8] for 
more applications. For the viscoelastic materials, the competition between the kinetic 
energy and the internal elastic energy through the special transport properties of their re- 
spective internal elastic variables makes the materials more untractable in understanding 
their behavior, since any distortion of microstructures, patterns or configurations in the 
dynamical flow will involve the deformation tensor. For classical simple fluids, the internal 
energy can be determined solely by the determinant of the deformation tensor; however, 
the internal energy of complex fluids carries all the information of the deformation tensor. 
The interaction between the microscopic elastic properties and the macroscopic fluid mo- 
tions leads to the rich and complicated rheological phenomena in viscoelastic fluids, and 
also causes formidable analytic and numerical challenges in mathematical analysis. The 
equations of the compressible viscoelastic fluids of Oldroyd type ([211 122]) in three spatial 
dimensions take the following form [3 \T7\ |2l] : 

Qt + div(£)u) = 0, (1.1a) 
{gu)t + div {gu u) - /xAu - + A)Vdivu + VP(^) = div(£»FF^), (1.1b) 
Ft-hu- VF = VuF, (1.1c) 
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where g stands for the density, u G the velocity, and F E M^^^ (the set of 3 x 3 
matrices) the deformation gradient. The viscosity coefficients /x > and A satisfy 

^ > 0, and 2^ + 3A > 0, 

which ensures that the operator — /iAu— + A)Vdivu is a strongly elliptic operator. The 
term P{q) represents the pressure for the barotropic case and is an increasing and convex 
smooth function of ^ > with > 0. The symbol ® denotes the Kronecker tensor 

product and means the transpose matrix of F. As usual we call equation (jl.lap the 
continuity equation. For system ()l.ip . the corresponding elastic energy is chosen to be the 
special form of the Hookean linear elasticity: 

W{¥) = \\Y\\ 

which, however, does not reduce the essential difficulties for analysis. The methods and 

results of this paper can be applied to more general cases. 

In this paper, we consider equations (jl.ip in the three-dimensional bounded domain 
C with sufficiently smooth boundary, subject to the following initial and boundary 

conditions: 

f(^,u,F)|i=o = (£'o(a::),Uo(a;),Fo(x)), a; G 17, ^^^^ 
\u|an = 0, 

and we are interested in the existence and uniqueness of strong solution to the initial- 
boundary value problem (jl.ip - (jl.2p near its equilibrium state. Here the equilibrium state 
of the system (jl.ip is defined as: is a positive constant (for simplicity, g = 1), u = 0, 
and F = I (the identity matrix in M^^^). We denote the perturbations of the density, the 
velocity, and the deformation gradient about the equilibrium by p, u, and E, respectively, 
that is, 

g=l+p, F = I + E. 
Then, system (jl.ip becomes equivalently 

pt + div(/3u) -|- divu = 0, (1.3a) 
((1 + p)u)t + div ((1 + p)u (g) u) - pAu - (/i + A)Vdivu + VP(1 + p) 

= div((l + /j)(/ + S)(/ + S)^), (1.3b) 
Et + u- VE = V\iE + Vm, (1.3c) 
with the initial and boundary conditions: 

{{p,\i,E)\t=o = {pq{x),\iq{x),Eq{x)) = {g^ix) - 1, uo(2;), Fo(x) - I), x € 

By a strong solution, we mean a triplet (p, u, E) satisfying (jl.3p almost everywhere with 



the initial and boundary conditions (II. 4p . in particular, (p, u, £■)(•, t) G W^''^ x {W'^''^) 
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do) in this paper. 

When the density is a constant, system (11. ip governs the homogeneous incompressible 
viscoelastic fluids, and there exist rich results in the literature for the global existence of 
classical solutions (namely in or other functional spaces with much higher regularity); 
see P El [H Uni [m [El [III [H [II] and the references therein. When the density is not a 
constant, the question related to existence becomes much more complicated. In [6j the 
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authors considered the global existence of classical solutions of small perturbation near its 
equilibrium for the compressible viscoelastic fluids in critical spaces (a functional space 
in which the system is scaling invariant), see also [23]. For the local existence of strong 
solutions with large initial data, see [5]. One of the main difficulties in proving the global 
existence for compressible viscoelastic fluids is the lacking of the dissipative estimate for 
the deformation gradient. To overcome this difficulty, the authors in [10] introduced an 
auxiliary function to obtain the dissipative estimate, while the authors in [12] directly deal 
with the quantities such as Au + divF. Those methods can provide them with some good 
estimates, partly because of their high regularity of (u, F). However, in this paper, we deal 
with the strong solution with much less regularity for (u, F) E W"^'"^ x W^'*^, q € (3,oo), 
hence those methods do not apply (more precisely, we can not apply the standard energy 
method in our situation). For our purpose, we first obtain a uniform estimate for the 
linearized momentum equation using the maximal regularities of the Stokes equations and 
parabolic equations, and then we find that a combination between the velocity and the 
convolution of the divergence of the deformation gradient with the fundamental solution 
of Laplace operator will develop some good dissipative estimates which are very useful for 
the global existence. 

When the initial-boundary value problem (jl.3p - (jl.4p is considered, several difficulties 
need to be addressed: 

(1) For the initial-boundary value problem, the boundary condition can not be pre- 
scribed as zero for the density or the deformation gradient, and hence taking the 
derivatives of equations up to arbitrary orders and integrating by parts can not be 
applied. To overcome this difficulty, we need to establish a uniform estimate for 
the linearized momentum equation with the help of the standard — L"? estimates 
for the Stokes equations and the variants for the heat euqation; 

(2) The functional framework now is the general L'^ space with q > 3. This setting 
obviously exclude the classical methods which are the key tools in [5l \T0[ [T3l [HI [23] . 
To handle this difficulty, the classical Poincare's inequality is used and a new 
conserved quantity for (jl.ip will be needed. 

One key observation of this work is that under the condition ()2.4p (see Section 2) on 
the curl of the deformation gradient initially, not only the curl of the deformation gradient 
at any positive time is a higher order term, but also is the sum of the gradient of density 
and the divergence of the deformation gradient (see (15. 4p ). which motivates us to use 
the — L'^ estimate of Stokes equations. Furthermore, the divergence of the difference 
between the deformation gradient and its transpose is also a higher order term, see (I5.24p . 

The viscoelastic fluid system p.ip can be regarded as a combination of the inhomo- 
geneous compressible Navier-Stokes equations with the source term div(pFF~'') and the 
equation (ll.lcj) . For the global existence of classical solutions with small perturbation 
near an equilibrium for the compressible Navier-Stokes equations, we refer the reader to 
[181 Uni EQ] and the references cited therein. We remark that, for the nonlinear invis- 
cid elastic systems, the existence of solutions was established by Sideris-Thomases in [26] 
under the null condition; see also [25] for a related discussion. 

The existence of global weak solutions with large initial data of (11. ip is still an out- 
standing open question. In this direction for the homogeneous incompressible viscoelastic 
fluids, when the contribution of the strain rate (symmetric part of Vu) in the constitutive 
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equation is neglected, Lions-Masmoudi in [T6] proved the global existence of weak solu- 
tions with large initial data for the Oldroyd model. Also Lin-Liu-Zhang in [13J proved the 
existence of global weak solutions with large initial data for the incompressible viscoelastic 
fluids when the velocity satisfies the Lipschitz condition. When dealing with the global 
existence of weak solutions of the viscoelastic fluid system (11. ip with large data, the rapid 
oscillation of the density and the non-compatibility between the quadratic form and the 
weak convergence are two of the major difficulties. 

The rest of the paper is organized as follows. In Section 2, we recall several intrinsic 
properties of the system (jl.ip and also show a new conserved quantity. In Section 3, we 
introduce the functional spaces and state our main results, including the local existence and 
uniqueness of the strong solution to the system (jl.3p - ()1.4p . as well as the global existence. 
In Section 4, we prove the uniform estimate for the linearized momentum equation. In 
Section 5, the main goal is to obtain a series of uniform estimates for the gradient of the 
density and the gradient of the deformation gradient. In Section 6, we establish some 
uniform in time a priori estimates on the dissipation of the deformation gradient and the 
density, and finally finish the proof of our main theorem. 



2. Some Intrinsic Properties of System (jl.lj) 

In this section, we recall some intrinsic properties of system (II. ip and introduce a new 
conserved quantity. First, for the system (jl.ip it was proved that (see Lemmas 6.1 and 
6.2 in Hu-Wang [6j, and also p3]). 

Proposition 2.1. Assume that (^j, u, F) is a solution of system (II. ip . Then the following 
identities 

£>det(F) = 1, (2.1) 
div(^>F^) = 0, (2.2) 

and 

FifcViFy = (2.3) 
hold for all time t > if they are satisfied initially. 

Note that (|2.3p can be interpreted in terms of the perturbation about the equilibrium 

as 

^^k-^ij ~ dxjEik = EijViEik — EikViEij. (2.4) 

Similarly, (j2.ip can be also interpreted in terms of the perturbation about the equilibrium 
as 

{1 + p)dei{I + E) = 1, 

which implies 

(1 + p){l + tTE + i [{tvEf - tT{E^)] + det E^ = 1 (2.5) 
since for any 3x3 matrix E, we have 

det(/ + E) = 1 + trE + ^ [(trE)'^ - tT{E'^)] + det E. 
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iiE = -p- piiE + + p) 



[iT{E'^) - [iTEf] -deiE 



(2.6) 



It is worthy to pointing out that by a similar argument, the constraint in the two- 
dimensional case between the perturbations of density and the deformation gradient has 
the form of 

iiE = -p- piiE - (1 + p) det E. 
Similar to the conservation of mass due to the continuity equation (II. lap , the quantity 

JQ, 

is also conserved, where Fj is the i-th column of the matrix F. 

Proposition 2.2. // QQ{¥Q)idx = and (1\\{qq¥~q) = 0, then for all time t > 0, 

/ Q¥i{x,t)dx = 0. (2.7) 

Proof. Indeed, from (jl.lap and (jl.lcp . we deduce that 

dt{QY) + u • V(£»F) = Vu(£»F) - ^iFdivu. (2.8) 

On the other hand, from the vector identity 

curl(A X B) = AdiwB - BdiyA + {B ■ V)A - {A ■ V)B 

for all vector-valued functions A, B, and div(^)F^) = 0, one has 

u • V(£»Fi) - Vu£>Fi + QYidiYU 

= u • V{Q¥i) — gFi ■ Vu + ^jFjdivu 
= -V X (u X gFi). 

Thus, we can rewrite (j2.8p as 

dtigFi) - V X (u X gFi) = 0. 
Integrating this identity over Q gives 

d f 

— J^gFidx = 0. (2.9) 
The proof is complete. □ 
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3. Main Results 

In this paper, the standard notations for Sobolev spaces W'^''^ and Besov spaces Bp^ wih 
be used, and the fohowing interpolation spaces will be needed: 

and 

Now we introduce the following functional spaces to which the solution and initial condition 
of the system (jl.3p will belong. Given 1 < p,q < oo and T > 0, we set Qt = x (0,T), 
and 

n;P.'?(0, T) := {u € W^'P{0, T- {Li{^)f) n LP(0, T; {W'^^i{^)f)] 
with the norm 



'^\\v\'P't{0,T) '■— l|Ut||LP(0,T;L9{n)) + l|u||LP{0,T;H^2,9(f^)) , 

as well as ^ 



with the norm 

IK/, 5) Hyp-" := 11/11 2(1-1) + 11/11 i_i + \\g\\w^,i{n)- 

VYp 1 p 

For simplicity of notations, we drop the superscripts p,q in W^''^ and Vq'"^, that is, we 
denote 

W := WP'^ Vo := Vjf'". 
For large initial data, the following local in time well-posedness can be obtained as in 

Theorem 3.1. Assume that ft is a bounded domain in M.^ with C^'^^ (P > Q) boundary 
and (uo, /Oo, -E-o) € Vb with p G [2,oo),q G (3,oo). There exists a positive constant Tq 
such that the initial-boundary value problem (jl.3p - (jl.4p has a unique strong solution on 
fl X {0,Tq), satisfying 

(u, p, E) G W(0, To) X {W^'P{Q, To; L^(l^)) n ^^'(O, To; W'''^m) '° . 

The argument for proving Theorem 13.11 is similar to that in [5j , thus we omit the proof 
here. 

Remark 3.1. An interesting case is the case q < p. Indeed, by the real interpolation 
method, we have 

on i\„ 2(i-i) 2(i-i) 

and 

p'^ CB,p^ =Zp ^ 



Then, if we replace the functional space Vq''^ in Theorem 13.11 by 
Theorem 13.11 is still valid. 
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Now our main result can be stated as follows. 

Theorem 3.2. Assume that Q, is a bounded domain in with C^'^^ (j3 > Q) boundary 
and {vLQ, po,Eq) G Vq with p G [2,oo),q € (3,oo). There exists a positive number R < 1 
such that if the initial data satisfies ()2.2p , (I2.3p , ()2.7p , and 

!|(uo,po,^o)llVo <R^ 
then the initial-boundary value problem (|1.3p - (|1.4p has a unique global strong solution 

iu,p,E) G >V(0,oo) X (Tyi'P(0,oo;L'?(!^))nLP(0,oo;T^i'«(!^)))^°, 

satisfying 

|u||w(0,oo) < R'l 
|p||L°°{0,oo;W^i.9{n)) < 



\E\ 



L°°(0,oo;H/l''J(n)) 



< R. 



The rest of this paper is devoted to the proof of Theorem 
tation, we will assume that 

\n\ = i, 



(3.1) 

For simplicity of presen- 



podx = 0, 



and 



/ 

Jn 



{1 + Po)iI + Eo)ijdx 



0, if j; 

1, if i= j. 



(3.2) 



(3.3) 



Note that from the continuity equation and (|2.9p in Section 2, the identities (|3.2p and 
p.3p hold also for all positive t > 0. 



4. Estimates for the Linearized Momentum Equation 
In this section, we consider the following linearized equation: 

dtf - pAf - (/i + A)Vdiv/ + Vh = g 

/(O) = /o, flan = 0, ^ ■ ^ 

where / G W(0,T), g G (LP(0, T; L'?(^])))^ h G LP{0,T; Li{n)), and /o G V^''^{n). For 
this equation, we have 

Lemma 4.1. Assume that Q is a bounded domain in with C"^^^ (jS > Q) boundary. 
Then there is a positive constant C , independent ofT, such that 

ll/llw(0,T) + l|V/l||LP(o,T;L9(n)) < C {\\g\\LP{0,T:L'i{^)) + II /o II \/(f '"(Q)) • (4-2) 

The strategy of proving Lemma l4. II is to decompose the function / into two parts: one 
is a solution to Stokes equations and the other is a solution to a parabolic equation. For 
this purpose, let us first recall the maximal regularity theorem for parabolic equations (cf. 
[20]). 
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Proposition 4.1. Assume that i7 is a bounded domain in with C^^^ (f3 > 0) boundary. 
Given 1 < p < oo, ujq G Vq''' and ip E LP{0,T; L'^{Q,)), the Cauchy problem 

J If - fiAu - (/i + A) Vdivtj = if, 
\uj{0) = (Jo, 

has a unique solution u € W{0,T), and 

ll'^llw(0,T) < Ci (\\^\\LP{0,T;Li{n)) + ll'^o||yP'«) , 

where Ci is independent of ojq and ip. In addition, there exists a positive constant cq 
independent of tp such that 

||a;|[w(o,T) > Co sup \\oj{t)\\yv,<i. 

t&(Q,T) 

The classical — estimates for Stokes equations can be stated as follows (see The- 
orem 3.2 in [3]). 

Proposition 4.2. Assume thatQ is a bounded domain inM^ with C'^'^^ (fi > boundary, 
and G L'^(IR+, L'^(il)) for 1 < q,s < oo. Then the initial-boundary value problem of the 
Stokes equations 

dfU — fiAu + Vvr = Ip, irdx = 0, 
divu = 0, 
u\oQ = u\t=o = 0, 

has a unique solution (u, vr) satisfying the following inequality for all T > 0: 

||(Au, V7r,Uj)|[^«(0,7..i,(Q)) < C'||V'||l«(0,T;L'J(C)) 

with C = C{q, s, Q, fi). 

Proof of Lemma \4.1\ To begin with, we first introduce two new functions /i and /2, which 
satisfy respectively 

'dtfi-fiAfi + Vh = g 
div/i = (4.3) 
Ui(0)=0' fi\dn = 0, 



and 



Notice that 



dtf2 - - (/i + A)Vdiv/2 = 
/2(0) = /o, f2\dn = 0- 



/ = /l + /2 

by the linearity of the equation (j4.ip . For Stokes system ()4.3p . using Proposition 14.21 there 
exists a solution fi to ()4.3p and it satisfies 

ll/l||w(0,T) + l|V/l||LP{0,T;L9(r2)) < ||5'||ip(o,T;L'J(n)) , 

where the positive constant C depends on p, q, Q and does not depend on T. On the other 
hand, by the maximal regularity theorem of parabolic equations, there exists a solution 
/2 (extended to zero outside 0.) to ()4.4p and /2 satisfies 

ll/2|lw(0,T) < C'll/o|IV|f'«> 
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where the positive constant C depends on /x, A,p, q and does not depend on T. Hence, we 
have 

l|/llw(0,T) + l|V/l||LP(o,r;L9(n)) < ll/l|lw(0,T) + ll/2||>V(0,T) + II V/l||LP(o,T;L9(n)) 

< C {\\g\\LP{Q,T;Li(^)) + II/olly(f'«) ; 

where the positive constant C depends on p, q, Cl and does not depend on T. 

□ 



5. Uniform Estimates 

Due to Theorem 13.11 for any given initial data, the local existence and uniqueness of 
strong solution can be established. In order to extend the local solution to a global one, 
we need to establish a series of uniform bounds. That is the main objective of this section. 

Throughout this section, we assume that over the time interval [0, T], for a given suffi- 
ciently small positive number R < 1, the following bounds hold 

u||w(o,T) < R; 

-E'||L°o(o,r;iyi.9(n)) < R- 

The main goal of this section is to obtain some uniform bounds on Vp and VE in 
LP{0,T-Li{n)) (see Corollary O below). 

To begin with, notice that since q > 3, ()5.ip will imply 

||p||L°°(f7) < C'll/'lltyi^9(n) < CR < 2' 
if R is sufficiently small. Similarly, one can assume 

ll-E-jj ll^^oo < CR < 1, for all i,j = 1,2,3. 

5.1. Dissipation of the gradient of the density. To prove Theorem 13.21 valid . we need 
first to establish a uniform estimate on the dissipation of the gradient of the density. 

Lemma 5.1. Under the same condition as Theorem \3.^ the solution {p,u,E) satisfies 

\\^P\\LP[0,T;Li(n)) < C{R^ + -R||V£'||lp(o,T;L9{Q))), (5.2) 

where C = C{p, q, fi, A, il). 

Proof. To begin with, we introduce the zero-th order differential operator 

V ■ ■ — A~^8 B 

which is defined as 

^^./ = (fl-^/) ' 

where J- denotes the Fourier transformation. Note that 

ll^ii/lli^ < C||/||l,, for l<g<(X) and 1 < i, j < 3. 
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Extending the perturbations of the density and the deformation gradient by outside the 
domain Q, applying the operator TZij to (]l.lbp . and using the equation (jl.lap . we obtain 

TZij ((1 + p)ut + (1 + p)u • Vu) - {2fi + A)a^,divu + ad^^p 

(5 3) 

= TZij {ptrE + (1 + p)Eikd,,Eik) + Mj + (P'(l) - P'{1 + p)) d^^p, 

where 

a = (l + P'(l)), 
and we used the fact that, from ()2.4p and ()2.6p . 

TZijdiYE = A'^d^^dx,d:„^Eik = A"^9^^.5a;^(9a;,£'jfc 

= A ^dxjdx^dx^Eii + TZjk {EiiViEik — EikViEu) 
= dxjt'rE + TZjk {EiiViEik - EikViEa) 

(5.4) 



-dx^p-dx^{piTE) + dx^ [{l + p) 

+ TZjk {Eii^iEik — EikViEii) 
-dxjP + Mj, 



^ [ti{E^) - itvEf] -detE ^ 



with 



i [ti{E^) - {tiEf] -detE^+ TZjk {EuViEik - EikViEu) . 



2 

We rewrite the equation (|5.3p as 

(JZiju)t — fiATlijU — {fi + X)VdivTZijU + ad^^p 

= -TZij {put + (1 + p)u ■ Vu) + n^j ipiiE + (1 + p)Eikdx^Eik) + Mj (5.5) 

+ {P'{l)-P'{l+p))dx,p. 

According to Lemma l4.H we have 

ll^iiU|lw{0,T) + a||Vp||LP(o,T;L9(n)) 

< c(\\TZijUo\\yP,i + WTZij {put + (1 + p)u • Vu) \\LP{0,T;Li{n)) 

+ \\TZij {ptrE + (1 + p)EikdxiEik) \\lp{o,T;Li{^)) 

+ \\Mj \\LPiO,T-Mm + II (^'(1) - ^"(1 + P)) 5x,P|L,(o,^^^,(^)) 

Next, we estimate the terms on the right hand side of the above inequahty. Indeed, 

\\TZij {put + (1 + p)u ■ Vu) \\LP{0,T;Li{n)) 

< \\put + (1 + p)u ■ Vu||ip(o,T;L'J{Q)) 

< \\p\\L°°\\ut\\LP{0,T;Li{n)) + (1 + II Pll oo) || u||loo (o,T;L'3(n)) II Vu|| £^p(0,T;L°° (H)) 

< (^(^IIpIIl"" ||ut||Lp(o,T;L'3(n)) + (1 + \\p\\od)\\^\\l°^ {0,T;Li {n))\\^\\LP {0,T;W'^-i{n)) ^ 



(5.6) 
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WUij {ptiE + (1 + p)EikdxiEik) \\LP{0,T;Lim 

< WptiE + (1 + p)Eikd^^Eik\\LP{o,T;Li{n)) 

< ||-E||L°°||p||LP(0,T;L9(n)) + (1 + ) ll-^'ll L°° II V£^||LP(o,T;L9(n)) 

< C {R^ + R\\VE\\L,^O^T;L'>m) 

\\'^j\\LPiO,T;Li{n)) < W^iEiiVlEik - ElkVlEii) ||lp{o,T;L9(Q)) 



+ 



- [tT{E'^) - [iiEf] -detE 



2 

< C lEik — EikV lEiiWi^pi^Q^j^.iq^p^)-^ 

+ P\\LP{<d,T;Li{n)){\\E\\\^ + ||-E||i=o) 

+ (1 + ||/0||L-)||Vi?||iP(0,T;L.(n))(ll^l|L- + II^IIL 

< c(i?||V^||iP(o 



LP{0,T;Li{fl)) 

(5.7) 



since R < 1; 



II (P'(l) - P'(l + p)) 5.,p|Lp(o,r;Lnf^)) - ll^-.^ILf{0,T;Lnf^)) ll^l 

< Ci?||V/9||iP(0,r;L9(O)), 

here we used the fact 

P'(l + /?) — -P'(l) = P"{z)p for some z between 1 and 1 + /o, 

and hence 

\\P\l + p)-P'{l)\\L^<ri\\p\\L^ 

with 

ry = sup |P"(2:)|. 

2 -- 2 

Substituting those estimates back into ()5.6|] . we obtain 

ll'^jju||w(0,T) + a||Vp||LP(o,T;L'?(n)) 

T;Li{n)) + R\\^P\\LP{0,T;Li{n)) 

Assuming that R is so small that CR < ^, we deduce from the above inequality that 

ll^iju||w(0,T) + ^l|Vp||Lp(o,r;L9(n)) < c(^R^ + R\\\7E\\LP{o,T;Liin))y (5.8) 
The proof is complete. □ 
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5.2. Dissipation of the deformation gradient. The main difficulty of the proof of 
Theorem 13.21 is to obtain estimates on the dissipation of the deformation gradient. This is 
partly because of the transport structure of equation (11.3cp . It is worthy of pointing out 
that it is extremely difficult to directly deduce the dissipation of the deformation gradient. 
Fortunately, for the viscoelastic fluids system (II. Sp . as we can see in O |9l [lOl [I2l [HI 
114] . some sort of combinations between the gradient of the velocity and the deformation 
gradient indeed induce some dissipation. To make this statement more precise, we rewrite 
the momentum equation ()1.3bp as, using (j2.2p 



and prove the following estimate: 

Lemma 5.2. Under the same condition as Theorem \3.S\ the solution {p,u,E) satisfies 



dfU — iiAu — divE 



-(l + p)(u- V)u- VP(l + p) 

+ pdivE + (1 + p)Ejkdx,Eik - pdtu, 



(5.9) 




(5.10) 




(5.11) 



Then, (15. 9p becomes 




(5.12) 



where 



(1 + p)(u • V)u - VP(1 + p)+ pdivE + (1 + p)Ejkd^^Eik - pdtU. 



Also, from (|1.3cp . we have 




C div(Vu + Vu^ - (u • V)^) 



(5.13) 



From (I5.12P and (I5.13p . we deduce, denoting Z = u — f^i, 

dtZ - pAZ -{p + A)VdivZ = n:=ni-'H2 



(5.14) 



where 



^2 = -/:(div(Vu + Vu^; - (u • V)^)). 



Equation (|5.14p with Proposition 14.11 implies that 




(5.15) 



< C{p,q) {R+ \\'H\\LP{o,T;Li{n))) ■ 
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Next, we estimate ||^i||LP(o,T;L9(n))! i = 1,2. Indeed, for T-Li, using (|5.2p . we have 

ll^l||LP(0,T;L9(n)) < (1 + l|p||L°°(QT))l|u|lL°°(0,T;L'J(n))||Vu||/,p(o,T;L°°(C)) 
+ «^I1Vp|1lp(0,T;L'3(C)) + \\p\\L^{QT)\\'^^\\LP{0,T;Li(n)) 
+ (1 + ll/'llL°°(QT))ll^lli°°(QT)ll^-^lliP(0,T;L9(n)) 

+ \\p\\L°°(QT)\\'^t^\\LP{0,T;Li{n)) (5.16) 
< C(^i? + -R|lu||ip(o,T;W2>'J(n)) + ^11 V-E'||LP(0,T;L'J(n))) 



< C 



since i? < 1 with 



For T-L2, we have 



,T;Li{n)) , 



K = sup P'(l + Z). 



|^2| < -|u| + -£{div{WuE - (u • V)^;)), 



and 

||Vu£^ - (u • V)£'||iP(0,T;L'?(Q)) < l|Vu||ip(o,T;L9(Q))l|-E'l|L°° 

+ l|V-E||LP(o,T;L'J{n)) 

< R^. 

Hence, one can estimate, by the standard estimates of eUiptic equations, 

II^2||lp(0,T;L9(Q)) < -||u||LP(o,T;L'?(n)) + -||Vui? - (u • V)i?||ip(o,T;L9(n)) 

< C{R + R^) < CR. 
Therefore, from ()5.16p and ()5.17p . we obtain 

\m\Lv{0,T;L'>m <c{r + R\\VE\\Lvio,T;L'>m) ■ (5-18) 

Inequahties (|5.15p and (|5.18p imply that 

ll^llLP(0,T;VK2.9(f7)) <c{r + R\\VE\\lv{0^T;Li(O.))^- (5.19) 

Hence, we have, from (|5.1ip 

||div£^|[iP(0,T;L9(f7)) < (||^||lp(0,T;W2.9{Q)) + W^W LP {0,T;W'^'i(Q.))) 
<C(i? + i?||Vi?||ip(o,T;L.(C)))). 

On the other hand, from the identity (|2.4p . we deduce that 

||curl -E'i||iP(o,r;L9(n)) < '^\\E\\L°°{Qj,)\\VE\\LP^Q^T;Li{n)) 

< C||-E'||L°°(0,T;iyl>'3(Q))||V£'||iP(o,T;L'J(0)) (5-21) 

< C-R||V£'||iP(o,T;L9(n))- 

Combining together ()5.20p . ()5.2ip and the inequahty 

||V/||l. < C(||div/||L. + llcurl/IUO for ah 1< r < oo, 



(5.20) 
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we obtain 

W^^W LP {0,T;Li{n)) < (^(^ + -R||V-E'||LP(o,T;L9(n))) ) 

and hence, by choosing CR < ^, one obtains (jS.lOp . 

The proof of Lemma 15.21 is complete. □ 

Notice that, from (j5.2p and ()5.10p . one can improve the estimates on the gradient of 
the density and the divergence of the deformation gradient. Indeed, we have 

Corollary 5.1. Under the same condition as Theorem \3.^ the solution {p,u,E) satisfies 

l|Vp||iP(o,T;M(c)) < CR^ (5.22) 

and 

I|V£'||lp(o,T;L9{Q)) < CR^, (5.23) 
where the positive constant C is independent ofT. 

Proof. First, note that the estimate ()5.22p is a direct consequence of Lemma 15.11 and 
Lemma 15.21 

To show ()5.23p . according to ()5.2ip and Lemma 15.21 we only need to verify 

||div£;||iP(0,T;L9(f7)) < CR^. 
To this end, we notice first that, since divcurl/ = for any vector- valued function /, thus 

(-A)"Mivcurldiv£; = 0. 
On the other hand, we also have, using the identity ()2.4p . 

(-A)-Mivcurldiv^ = {-A)-^dw{d^^d^^Eij - d^^d^^Ekj) 
= {-A)-'^div{dx^dr,^Eik - d^^d^^Eki) 

+ {-A)-^diYd^^{EijViEik - EikViEij - EijViEk, + EiiViE^j) 
= divE - divE^ + Af 

with 

M = {-A)-^dwd^^{EijViEik - EikWiE.j - EijViEki + EuViEkj). 
Hence, we have 

div^; = div^"^ - A^. (5.24) 

The property (|2.2p implies that 

divE^^ = -V/9 - dw{pE'^), 
which together with (j5.22p yield 

l|div£'^|liP(o,T;L9(n)) = W^P + (i^'^{pE~^)\\LP{o,T;Li(n)) 

< \Np\\LP{0,T;Li(n)) + l|V£'||iP(o,T;L'J(n)) 
+ \\E\\L^\\^P\\LP{0,T;Li(n)) 

< CR^. 



(5.25) 
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For TV, we have, using Lemma 15.21 

(5.26) 



LP{0,T;L1{Q)) < 4||£'j|Loo || Vii^|| j^p(o,T;L9(n)) 



< CR^ 

Combining (^^M), (fOSj) . and ([Oe]) . we have 

\\divE\\Lp(^Q^j-,Li(n)) < CR^. 
The proof is complete. □ 

6. Proof of Theorem 13.21 

The goal of this section is to prove Theorem 13.21 To this end, for a given R > small 
enough, we define 

Tma.x '■= sup |t > : there exists a solution (p, u, E) to (|1.3p with ||u|[w(o,T) < R aiid 

^^^{\\p\\L°°(0,T;W^'i{n))7 ll^llL°={0,T;H/i'9(n))} < 

We first show that in the interval [0, Tmax]) we have better estimates on the velocity. 

Lemma 6.1. For any T G [0,Tmax], we have 

l|u|lw(o,T) < CR^ 
for some positive constant C independent ofT. 

Proof. In fact, by Proposition 14.11 and equation ()1.3bp . we have 

||u||h;(o,T) < C{p,q,n)(^\\uo\\yP,i + \\pdtu\\Lp(^Q^T;Li{n)) 

+ 11(1 + Vu||ip(o,T;L9(n)) 

+ l|VP(l + p)||LP(o,T;Lnf^)) (6-1) 
+ 11(1 +P)divS||Lp(o,T;L9(n)) 

+ 11(1 + P)EjkdxjEik\\LP(0,T;Li{n))^ ■ 

For those terms on the right hand side of (|6.ip , in view of Lemma 15.21 and Corollary 15.11 
we have the following estimates: 

||/?(9tU||iP(0,T;L9{Q)) < \\p\\L^\\dt'^\\LP{0,T;Li{n)) < CR^; 
11(1 + P)u • Vu||2,P(o,T;L9(n)) < (1 + IIpI1l°°)I|u||l°°{0,T;L9(Q))||Vu|1^p(0,T;L°°) 

< C||u||^oo(o,r;L9(n)) l|u||LP(o,T;iy2>'j(n)) 



< CR^; 

|VP(1 + p)\\LP{0,T;Li{n)) < l^\\'^P\\LP(0,T;Li{n)) < CR^] 



11(1 + /))div£;||ip(o,T;L'J{Q)) < (1 + l|p||L°°)||div£;||ip(o,T;L9{Q)) < CR^] 

(1 + p)£;jfc53,j.£'jfc||iP(o,r;L9(n)) < (1 + l|V£^|lLP(o,T;L9(n)) 

< CR^. 
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Substituting those estimates back into (j6.ip . we have 

||u||w(0,T) <C(||uo||^P.. + i?2) <Ci?2, 

provided that we choose the initial data to satisfy 

II^oIIvq''' — CR^. 

The proof is complete. 

□ 

According to Lemma |6.H we deduce from the above estimate that 

Next, we improve the estimates on the density and the deformation gradient. Due to 
the similarity of equations for the density and the deformation gradient, the arguments 
for those two unknowns are similar, and hence for the clarity of the presentation, we will 
only focus on the argument for the improved estimates for the density. To this end, we 
introduce a new variable: 

o" := V In g, 

and we have 

Lemma 6.2. Function a satisfies 

dta + V(u • a) = -Vdivu, (6.2) 

in the sense of distributions. Moreover, the norm \\cr{t)\\Lq(^Q) is continuous in time. 

Proof. We follow the argument in [20' (Section 9.8) by denoting as = S^o', where is the 
standard mollifier in the spatial variables. Then, we have 

dtae + V(u • as) = -Vdivu^ + 7^£, 

with 

Ke = V(u • ae) - SeV{u • a) 

= (u • Vae - Se{u ■ Va)) + (a.Vu - S,{a ■ Vu)) (6.3) 

=: nl + nl 

Since a G L°°(0, T; L''(0)) and u G LP(0, T; tyi'°°(J7)), we deduce from Lemma 6.7 in [20j 
(cf. Lemma 2.3 in [15]) that 

nl^O as e^O. 

Moreover, 
and 

Se{(T -Vu) ^ a -Vu in L^(0, T; L'?(J])) 
since a ■ Vu G L*'(0, T; L'?(il)). Thus, we have 

7^2^0 in LP{0,T;L''{n)). 
Then, taking the hmit as e ^ in ()6.3p . we get (16. 2p . 
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Multiplying (j6.2p by |cr|'^~^cj, and integrating over we get 



1 


^ II ll« 






Q 






/ 

Jn 



f (— Vdivu • cr|cr|'^ ^ — 5jUfc(Tj(Tfc[(T['' ^ divu|cj|'^)(ia 

Jn Q 



< \\u\\w2.,(n)\W\\%ln) + l|Vu||Loo||cj||f^, + ^||divu||L.c 

< (7||u||ty2,g||(T|||~^(l + ||0-||l9). 

q-1 



Dividing the above inequality by ||o"||^, , we obtain 



d .. „ 



Since a G L°°(0, T; L<?(17)), we have 



di 



< C\\u\\]y2,g{l + \\(t\\li] 



|ct||l, eLP{0,T). 



Thus, ||o"||l9 € C(0, r). The proof of Lemma [6^ is complete. 

With the aid of Lemma 16.21 we can improve the estimates on the density. 
Lemma 6.3. For any T G [0,Tmax], the density g satisfies 



liv^ll 



L°°{0,T;Li{n)) 



< cm 



where the positive constant C is independent ofT. 
Proof. We only need to prove 

3 

W\\L°°{0,T;Li{n)) < CR2. 

Indeed, multiplying (j6.2p by |(T|''~^cr and integrating over Q, yields 
1 d „ 



/ Vdivu • o-|o-|«"2 - / V(u • cr) • crlcrl^-^dx. 
Jn Jn 



Note that 



and 



where 



Vdivu ■ a\a 



19-2 



< \\^-^\\LHn)M%lny 



/ V{u-a)-a\a\^ '^dx = / dx^Ukakcrj\a\^ "^dx + I 
Jn Jn 

I =■■ d^i^dx^dx^{ln g)dx. In Q\a\'i^'^dx 
Jn 



aWWY'-'dx 



1 

q Jn 



^kdx^.{\o'\'')dx = — / |cj|''divudx 



QJn 
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Hence, 

V(u • a) ■ a\a\'^~'^dx 



< C||Vu||Loo||(7||^q < C||u||yi/2,q(f^)||(7||^, 



Substituting those estimates back into (j6.4|) . one obtains 

- C'lkllL9^l|u|lH/2.'j{Q) (1 + ■ 

Multiplying the above inequality by integrating over [0, T], we have, using 

Corrollarv 15.11 and Lemma |6.1[ 

-\WiT)\\Li < -Ikolliq + C'lkllip(o,T;L9(n))ll"llif{0,T;iy2,'J(n)) (l + Ik ||loo (o,T;L'J(Q))) 

< i||ao||i, + CR'P (1 + ML^io,T;LHn))) ^^'^^ 

< CR^P (1 + \W\\L^{0,T;Liin))) ■ 

Thanks to Lemma 16.21 the continuity of t i— > ||o'(t)||2,<7(n) implies that 



in some maximal interval (0, T(i?)) C (0,Tmax)- If T{R) < Tmax, then we have 

lk(r(i?))||i,(o) = i?t. 

On the other hand, (j6.5p implies 

R-2 = \HT{RmL,^n) < CvR^ (l + ||a||ioo(o,r;Lnn))) 
< CvR^{l + Rl). 
This is impossible, if we assume that R is so small that 

CpR^{1 + R^) < 1. 

Thus, this implies that 

3 

lFllL°°(0,Tmax;L«(O)) < R^ ■ 
The proof is complete. □ 

Similarly, we have 

3 

!|V£;|lLoo(o,T„^^^;L9(n)) < CR2, 

and hence 

3 

max{||Vp||Loo(o,r^^^;L.(n)), || V^||loo(o,t^^,;L9(q))} < CR2. (6.6) 

On the other hand, Proposition l2.2( combined with the initial conditions (|3.2p and (13. 3[) . 
imply that 

/ pdx = 0, and / (1 + p)Eij = 0, for aU i,j = 1, 2, 3. (6.7) 
Jn Jn 

Poincare's inequality, combined with (I6.6P and (I6.7p . yields 

3 

max{||p||ioo(o,T^^^;L,(n)), ||-E^||L°o(o,T^,,;L9(n))} < CR^, 
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and thus 

3 

^^^{\\p\\l°° (o,Tm^^;W^'i(n)) ^ ll-£^llL°°(o,rn,ax;VKi''?{n))} < CR2. (6.8) 
Now we are in a position to finish the proof of Theorem 13. 2i 

Proof of Theorem \3.S\ Suppose that Tmax < oo. Let T„ Tmax be arbitrary. Then 
necessarily 

l|u||w(o,T„) R as n ^ oo. 
Indeed, if this was not the case, then 

sup||u||w(o,T„) < R 

n 

and the function u, defined on (0,Tmax)i would be a solution to the equation on 
(0,rjnax)- Taking u(rinax) as a new initial condition, from Lemma WA\ we get u(rmax) G 
V^''^. And then, according to Theorem 13.11 there exists a Tq such that there exists a 
unique solution to (jl.ip on (Tmaxi Tmax + Tq). This means we can extend the solution we 
obtained outside the interval (0,riiiax) to (0,Tmax + To) with ||u||w(o,Tmax+To) < R and 

niax{||p||ioo(o,Tm^^+To;Vl/i.9(Q)), ||-E'||L°°{0,r„,ax+To;Vl/i.9(n))} < R 

for some Tq > 0. This is a contradiction with the maximality of Tmax- Thus, it is impossible 
to have 

^max ^ CX). 

Hence, the solution is defined for all positive time. 
Finally, from the estimate ()6.8p . one has 

l|p||L°°(0,oo;l^l.9(f7)) < R 

and 

\\E\\l°°(Q,oo;W^'1{Q.)) < R, ■ 

The proof is complete. □ 
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